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ABSTRACT 


The atmospheric effect on the upward radiance of sunlight 
scattered from the earth-atmosphere system is strongly influenced by 
the contrasts between fields and their sizes. For a given atmospheric 
turbidity, the atmospheric effect on classification of surface 
features is much stronger for nonuni form surfaces than for uniform 
surfaces. Therefore, the classification accuracy of agricultural 
fields and urban areas is dependent not only on the optical 
characteristics of the atmosphere, but also on the size of the surface 
elements to be classified and their contrasts. Atmospheric 
corrections that do not account for the nonuniformity of the surface 
have only a slight effect on the classification accuracy, in many 
cases, while in other cases the classification accuracy decreases. Ii 
this paper, the radiances above finite fields are computed to simulab 
radiances measured by a satellite. A simulation case including 11 
agricultural fields and four natural fields (water, soil, savanah, an 
forest) is used to test the effect of the size of the background 
reflectance and the optical thickness of the atmosphere on the 
classification accuracy. It is concluded that new atmospheric 
correction methods, which take into account the finite size of the 
fields, have to be developed in order to improve significantly the 
classification accuracy. 
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THE EFFECT OF FINITE FIELD SIZE ON 
CLASSIFICATION AND ATMOSPHERIC CORRECTION 

1. IMTRCXXJCTION 


Efforts have been increasing during the last decade towards 
applying satellite inuiges of the earth to classify and analyze surface 
features. ITiese programs have been successful enough to enccureige 
development of improved satellite sensors such as the Itiematic Mapper 
and the ^^Itispectral Resource Sampler (Schnetzler and Thompson, 

1980). The earth's atmosphere has a significant beari.g on the 
interpretation of the images. The initial studies of scattering and 
absorption effects in cloudless atmospheres assuned that the surface 
was uniform (Herman and Browning, 1975; Otterman and Fraser, 1976; 
Fraser et al., 1977; and Horvath et al., 1980). Attempts were made to 
develop corrections, again assiining that the surface was uniform 
(Turner, 1972; Fraser, 1974; Potter, 1976; Turner, 1978; and Slater, 
1980). These methods depended on values of the aeroso) optical 
thickness, which has been measured from satellites with success over 
large bodies of water (Griggs, 1973, 1974, and 1979; Mekler et al., 
1977), surfaces with discontinuities in reflectance (Kaufman and 
Joseph, 1981), but unsuccessfully over land. Pearce (1977) showed, 
however, that nonuniform surface reflectance patterns outside of an 
instantaneous- f ield-of-view (IFOV) modify the effect of the 
atmosphere on the radiance, (see also Kawata et al., 1978; Haba et 
al., 1979; Otterman and Fraser, 1979; Mekler and Kaufman, 1980; and 
Dave, 1980). Schowengerdt and Slater (1979) recently concluded that 
the influence of neighboring pixels causes only a secondary 
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stjnospheric effect. In the works of KswAts et al. (1978), Haba et al. 
(1979), and Kaufnwin and Joseph (1981), atmospheric corrections were 
applied, taking into account the finite size of the fields. Prom the 
Monte Carlo calculations of Paarce (1977), or from the approximate 
calculations of Mekler and Kaufman (1980), it is clear that the 
atmospheric effect is much stronger in the case of finite fields, than 
in the case of infinite fields (for the same atmospheric conditions) 
and thus, it is definitely not a secondary effect but a primary 
effect, v^erever significant contrasts in surface reflectance exist. 

'Ihe purpose of the present work is to show that classification 
errors are usually much bigger in regions consisting of finite fields 
cf dimensions less than a few kilometers compared to errors vhere the 
fields are larger. In order to make accurate classification, the 
measured radiance should be corrected for the a'oi'ospheric effect. 
Atmospheric corrections that do not take into account the finite size 
of the fields, reduce the r^ldiance errors slightly in most of the 
cases, in some of the cases the correction increases the error, and in 
the remaining few cases, the correction is satisfactory. 

The radiances given here are computed only for the nadir, for 
nunabsorbing atmospheres, for the visible through near infrared 
spectral band (400 nm < X < 1100 nm) , aixJ for a total normal optical 
thickness of < 1.0. Each portion of the earth's surface is assuned 
to reflect light isotropically according to Lambert's law. The term 
"field" here shall apply to either the test or the reference area, 
which is of uniform reflectance (e.g. a corn field a lake or an urban 
area). Ihe background is the region surrounding the field. 
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topics that follow ars arranged such that the basic surface 
reflectance characteristics of finite fields are reviewed in section 
2. Ihen a discussion of tSie atjnospheric effects associated with 
finite fields follows in section 3. Ihese effects cause errors in the 
derived iSdiances and attempts to r,ake corrections without accounting 
for the surface nonuniformity are discussed in section 4. Section 5 
contains a study of the classification errors in nonuniform regions. 
Ihe conclusiontN are given in section 6. Appendix 1 gives a list of 
the notations. Ihe equations used in these calculations are explained 
in the other appendices. 

2. SURFACE AND ATMOSPHERIC CHARACTERISTICS 

A. Surface reflectance 

Statistics of the size of agricultural fields and vMter bodies 
were recently reported by Pitts and Badhwar (1980). They shoved that 
the peak of the distribution of field size (including all crops) 
occurs for size of 6 acres (corresponding to a square of edge ler>gth 
of 155 m) , while 50 percent of the crop area is in fields smaller tnan 
130 acres (corresponding to a square of edge length of 720 m) . 
Therefore, a representative field-edge-length is of the order of one 
tenth to a few kilometers. Fbr such surface cover, it is necessary to 
use atmospiieric models which account for the finite size of the 
irtderlying fields. Figure 1 shows an example of the variations in the 
surface reflectance, taken by a Thematic Mapper simulator (Labovitz et 
al., 1980), and kindly given to us by D. L. Tbll. The image of Fig. 1 
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Figure 1 . An image of itie eastern fringes of Denver, Coiorado. acquired by • Thematic Mapper simuiator shows the hi^ rtonuniformity of the surface 
The area pee is 4 x 8km. The wa«e<e''gth is 450-520am, artd the terre is 1300 MOT (1900 GMT). 
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VMS acquired on the plains at the eastern fringe of Denver » Colorado. 

flight took place on June 20, 1979, starting at 1:00 P.M. MDT 
(1900 GMT). Ihe InuK^e is taken in the 450-S20 rvn band. Ihe 
underlying surface is an urban and rural area. Ihe image size is 4 x 
8 km; the size of the fields range betvMen a fev^ tenths to a fev^ 
hundred metets. 

Figure 2 shows the probability of reflectivity of the data given 
in Fig. 1, in four bands (450-520 nm, 520-600 nm, 630-690 nm, and 
760-900 rm) . The range of the reflectance extends from 0.02 to 0.34. 
Ihe reflectance was calculated by using laboratory calibration of the 
instrunent (Richard, 1979 and 1981). 

B. Atmospheric characteristics 

The atmospheric optical characteristics depend strongly on the 
aerosol optical thickness. Flowers et al. (1969) reported the results 
of a five year study of a . nospheric optical thickness measurements in 
the Uhited States. They found for a wavelength of 500 nm that in most 
of the cases the aerosol optical thickness for rural regions was in 
the range 0.1 £ £ 0.2, for suburban regions in the range 0.1 < 

0.3, and for urban regions 0.4 < t^. A recent study (Pearson et al., 
1981) analyzed turbidity data for 6 years from a non-urban area in 
North Carolina and shovMd an average optical thickness of ■ 0.3. 

'Hie simmer aver£xge was about ■ 0.6. These results show the h^.gh 
turbidities in the eastern U.S. 
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PROBABILITY RJNCTION 
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APPARENT REFLECTANCE 

Figurt 2. Tht pfobcbilitY of rtftoctivtlv of tfit data from tht imaga on Figura 1 tn fojr barKfi. 



NONUNIFORM SURFACE REFLECTANCE 


Figurt 3. Schamatic diagram of tht oontributiom to tht upward radiarrct. 
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3. THE EFFECT C3P THE FINITE SIZE 


In ord«r to discuss th« offset of nottunif^nn surfacss on ths 
radlancs msasursd at a satsllits, it is convsnisnt to uss a normalizsd 
radiance (X) , equivalent to the apparent reflectance, which is related 
to absolute radiance (I*) by 

I - "IV(^„F^), 

o o 

vi^ere • cos ® ® is the zenith anale of the sun; and F^ is the 

o o o o 

solar eoectral irradiance above the atmosphere on a unit area 
perpendicular to the beam. Ihe radiance of light reflected from the 
earth and iis atmosphere <s composed of three components (see Fig. 3): 


I 


^0 


( 1 ) 


where, 

Iq - radiance of light scattered from the direct sun beam into 
the detector's f ield-of-view without being scattered by the 
surface. Ihis component is independent of surface 
reflectance ar.d causes a loss of contrast on the satellite 
image. 

Ig - radiance of light reflected by the observed field and 

directly trananitted through the atmosphere. This component 
is the attenuated signal, the measurement of which gives 
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ramot* Mnalnq information about tha surfaca. 

- radl^nca of light raflactad from tha flald and Its 

background and than la acattarad by the atmoaphiia to tha 
datactor. Ihla componant modlflaa tha apparant apactral 
algnatura of tha obaarvad flald and raducaa tha contraat 
between It and tha adjacant rag Ion. 

'Hie diffaranca between the atmospheric effect for uniform and 
nonuniform surface reflectance can ba explained by means of these 
tnree components. Since tha term is independent of surface 
reflectance, it is the same for the uniform and nonuniform surfaces. 
Ihe term Ig is dependent mainly on the reflectivity of the observed 
field and therefore is almost ths same for uniform and nonuniform 
surfaces. ITie expression for Ig can be deduced from Chandrasekhar 
(1960, Bq. 201), where the explicit dependence of the upward radiance 
on the surface reflectance is given. Ihe resultant expression is: 

Y -T 
. o e o 

's • I - U *s (2> 

a 


where ^ is the total (diffuse and direct) trananission of the 
o 

sunlight to the surface, is the total atmospheric optical thickness 
(aerosols and molecules) Ag is the reflectivity of the field, and 
(1 - ) accounts for the contribution of multiple scattering of 

fl 

light between the surface and the atmosphere. The reflectance A^ is 
an avRr^lge of the reflectance of the field and the reflectance of its 
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background. Sine* th« tarm (lA ) is usually small, £A « 1 (for t 

A AO 

■ X ■ 550 ran, aixi A. ■ 0.2; tA. ■ 0.025). Iq dapsrvia mainly on 

AS A 

Ag rathar than <>n tha siza of tha fiald and tha raflactanca of its 
background. 


l^a diffararvea batvman tha atmospharic affact for uniform and 
nonuni form surfacas is exprassad by tha diffusa radianca, which is 
also daducad from Chandrasakhar (1960) i 




1 - 


(3) 


where D is tha diffusa traiisnission through tha atmosphera, rathar 
than diract transmission from tha surftca to tha datactoi; A^ is a 
weightad avarage (diffarant from of tha surfaca raflactanca of tha 
antira surfaca, including the ot'.arvad fiald and its surrounding. For 
a uniform surfaca, boch the fiald and its background are assuned to 
hm/e tha same reflectance, while in tha nonuniform case, they have a 
diffarant reflectance, resulting in different effects on the measured 
radiances. Itie values of A^ and A^ can In principle be computed by 
solving the equation of radiative transfer for tha nonuniform surfaca. 
In tha present work, the values are derived from Monte Carlo 
computations. 


Ibe term adjacency effect specifies tha influence of nonuniform 
areas on the radiance and is defined by the difference between 
radiances above uniform (I^) and nonuniform (I^) surface: 
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61^ - l" (Aj, Ag, T) - (A^, Aj, T) 


(4) 


Here, A^ and Ag are the reflectances of the observed field and its 
background, respectively. T is the atmospheric tu^Mdity. Hence, the 
adjacency effect is zero for a uniform surface, negative for a darker 
background, and postivie for a brighter backgrourK). 

The effect of a change in the atmospheric haziness on the 
detected radiance for uniform and nonuniform surfaces can be seen in 
Fig. 4. This figure gives a schematic representation of the changes 
in the three components of the radiance, for three different cases of 
background reflectance relative to the field reflectance: (a) 

background much brighter than the field, (b) uniform surface (s«ne 
brightness of the field and the background) , and (c) darker 
background. The effect of the atmosphere on two of the components, Iq 
and Ig, is the same in all of the cases. The diffuse radiance Ig 
decreases as the background reflectance decreases. The bright 
background in Fig. 4a causes a strong increase in the total radiance 
(I) with the optical thickness. The increase is moderate in Fig. lb, 
and the dark background in Fig. 4c causes a decrease in the radiance 
with the optical thickness. We conclude from this discussion that for 
the s-»ne atmosphere the radiance may increase or decrease, depending 
on the relative brightness of the field and its background. This 
phenomena can not be predicted by an atmospiieric model that is based 
on the assumption of uniform surface reflectance. 

Figure 5 shows the results of Pearce's (1977) Monte Carlo 
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OPTICAL THICKNESS 



OPTICAL THICKNESS 



OPTICAL THICKNESS 


Figuf# 4 Schcmatir dtagtam of th* atmui()hcric affact on tha 
upa«a((1 ladianca Both tba ladianca atid tha optical thicknati 
ara fl««an in arbitiary unitf *a| tha background ii much 
brightai than tha claMifiad araa; |b| uni*orm tuifaca (tama 
btighttmi). and (c) tha background ii much darkar. 


II 






computAtlona of the nwlir rAdl^nre Jibove the center of a tiquAre 
surroirtded by a different bAcKqro'TH (the IPOV is 30 m) , "Hie 
reflectAnce of the !it.|itAre ia ■ 0.4 for the (.ipper three curvea, 
vi4ille the bAckqiourvI ia of reflectance ■ 0.2. Fbr the lovwr three 
curvea, the aqiiare la of ref lei'tanc'e 0.2, wiiile the backqrourvl ia of 
reflectance A^. • 0.4. In thia fiqure, the railiarK'e ia ahown as a 
fitKtion of the edqe lemith of the a«.|uare. Tbe atimrsphere ia 
nonabaorbinq ami ia imxleleil with three different aerosol optical 
thicknesses (0.0, 0.213, ami 0.638), correapondlnq to turbidities (0, 
N, and IN) , respectively. 13ie aerosol size distribution is qiven by 
l^arce (1977). 'Hie imlex of reflectam'e ia m • 1.5 ♦ O.Oi. Tl>e 
asymptotic values of the rmliam'es are qiven on the riqht side for 
each case. It is seen that a chanqe in the atjnospheric turbidity from 
0 to 3N causes an increasie in the ratllance by 4 percent for an 
infinite m)uare (uniform surface) , Uiile for a finite s«.)uare (of 
reflectance ■ 0.4 on backqroumi of reflectarice ■ 0.2) of miqe 
lenqth of 0.3 km, the effect is to redui'e the railiance 15 percent wtien 
the aerosol optical tlrickness is increaseil fnmi (0.0) to (0.638). 

Tdris comparison emphasizes that the atmi'*spher ic effect is different 
for finite fields than for infinite fields. Fiqure 6 shows similar 
results for a sc^uare and backqrourvi of reflectancos 0.0 ami 0.6, 
respectively. This fiqure and Fiq. 5 are useil for the computations in 
this paper. 

In order to increase our confidence in the results of Fiqs. 5 and 

6, the results of the ftonto Carlo computations are compared to exact 

> 

miner ical results based on codes developeil by Dave (1970) for a 


12 


RADIANCE (I7Fo) 


EDGE LENGTH (KM) 



Figure 5. Upward radiance at The rtadir above the center of a touare as furwtion of the edge length of the square (afear Pearce. 1977). The atmoaphanc 
turtwlity it given for each curve The wavetength it 550 nm. The square rrflactnntv (fwld albedo) (Aj) arsd tfte background refteennee (A,) 
are ^ven for each group of curvet. The aav^nprotic radiartoe arsd *he radiance for infmitctimai aouare are compared with exact 

computations according to Dave's code Idathed arrows). 


irtiform surface, l^e comparison wns performed for the asymptotic 
radiances and for an infinitesimal square on a larqe background. In 
the later case, the results of Dave Mere used w.'th the aid of 
equations given in Appendix 2. Most of Dave's calculations ^lgree with 
the Monte Carlo results of Pearce (1977) with a difference of leas 
tlwin 2 percent. TTie Monte Carlo results M»re also compared with an 
analytical approximate solution by Mekler and Kaufman (1980), 
resulting in similar good aijreement. Itterefore, the results of Pigs. 

5 and 6, and similar data given by Pearce (1977) a^g «ised in the 
following analyses with confidence in their aoruracy. In order to use 
these computational results for different values of surface 
reflectance, the radiance transformations, explained in Appendix 3, 
are applied. 


4. THE CHANGE IN RADIANCE 


In this section we shall analyze the changes in the radiance used 
for classif ication due to the atmospheric effect in the case of finite 
fields. The classification is performed by comparing the radiance of 
the test area (I^) with the radiance above the reference area (I^.)* 

We shall use the following notation for the radiance 1^ (Aj , Ag^ , T^, 
Xj), where i can be t or r; A^ is the reflectivity of a square (the 

is the reflectivity of the background, 
is the atmospheric turbidity, and Xj is the edge length of the area. 
TTius, the radiance above the center of the reference area in Fig. 7, 
for example, is I^ (0.2, 0.1, 1, 0.29). The relative change in 
radiance is calculated by; 


reference or test area) , Ag^ 
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ASYMPTOTIC DAVE 
RADIANCE RESULTS 

1.6 

1.0 


.0 

.6 


3N— — — 3N 




I (A , A , T , X ) 
r r Br r r 


( 5 ) 


I 


(A 

e 




Where I (A^) is the radiance vi^ich would be measured above an area of 
reflectivity A^, v^en no atmosphere is present. Me shall always use 
the same reflectivity for the reference area and for the test area 
(A^ ■ A^) in all following examples; thus without an atmospheric 
effect, we would expect AI ■ 0. 


ITie effect of size and backgrourvi reflectance on the error in 
radiance is examined in Fig, 7. The surface and atmospheric 
characteristics of the reference area are the same in Figs. 7 and 8: 

A^ ■ 0.2, but its background reflectance A^^ ■ 0.1, ■ 0.29 km, and 

■ 1. The atmospheric turbidity is the same for the test area anti 
the reference areas in Fig. 7. Therefore, without an effect due to 
the finite sizes of the fields, the difference in the radiance vould 
be zero. It is seen that for the same background (A^^ • 0.1), the 
difference in the radiance for a small test area (X^ < 2 km), is small 

(|^| £ 0.025) but reaches - 0.075) for an infinite test area, 
e e 

For high reflectance of the background (Ag^ ■ 0.4), the difference 
reaches (0.25), while for a black background (Ag^ ■ 0), the difference 
reaches (-0.10). An atmospheric correction, which ignores the finite 
size of the fields, even if taking into account exactly all 
atmospheric optical characteristics, would not perform any correction 
on the data of Fig. 7, since the atmospheric optical thickness is the 
same for the test area and the reference areas. 
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EDGE LENGTH (KM) 

Figur* 7. Th« rttativt difftrtnc* in radiancn (lolid linat) of tba ta«t and raftrtnca areat (I, - at a 

function of tht Itngdt of tha tatt iquara iida (abacitM) and of tha raflactanca of tha background (Ag,) 
tunounding tha tatt tquara Tha rafararKa and tait araa hava tha tanta raflactarKa of A, ' A, ■ U.2, but tha 
raflactarKa of tha background to tha rafaranca araa it Ag, ■ 0.1. Tha rafaranca araa't adga langth it 
X ■ 0 29 km. Tha atmoipharic turbidity it tha tama for tha rafarartca and tatt araat 'T, ■ T, ■ 1). Tha 
dottad linat giva tha tolarancat for tha diffaranca in radiartca for danification purpoiat. 
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'IV«) tolerance differences in radiances are introduced in Pig. 1 , 
as well as in all other figures giving results for the difference in 
radiance. These two tolerances are representative for the maximum 
difference in the radiance for vhich classification is still possible 
(assuning there are no other sources of errors). In this case, the 
two tolerance values of ^ 0.05 and ^ 0.10 are chosen as the 
representative class separation in the classification examples given 
in section 5 of this paper. With such tolerance limits, correct 
classifications can be made when the b^lckground reflectance is 0.1 but 
not for 0.4 (Pig. 7) . 

The relative difference of radiances, vihen the reference and test 
fields are subject to different atmospheric turbidities is shown in 
Fig. 8. The turbidities for the test and reference areas are ■ 3 
and ■ 1. When the backgrounds to the test and reference areas have 
the same reflectances (Ag^ ■ Ag^ ■ 0.1), the relative difference in 
r^Kiiances ranges between -0.01 to 0.06 for fields smaller than 2 km. 
The magnitude of the radiance difference jumps to the 0.3 to 0.6 range 
for a large difference in the b^K:kground reflectances (Ag, ■ 0.4). 

The atmospheric corrections used in Fig. 8 for the uniform 
surface were derived from Fig. 9. There the upward radiance at nadir 
above a uniform s*jrface of reflectance (A) was plotted for the three 
atmospheric turbidities used in Figs. 5-8. The results shown on Fig. 

9 were calculated by Dave's routine for the same atmospheric 
characteristics that were used for the Monte Carlo computations. For 
"measured" radiance and turbidities appearing in Fig. 8, the 


18 



i/r7 



0.0 1.0 2.0 00 


EDGE LENGTH (KM) 

Figurt 8. Sam* ai in Figur* 7, but for diffarant atmoaphanc turbidity of tha tatt tnti rafararx:* araat 
T, ■ 1. T, ■ 3). Tha »oiid linai giv* th* ralativ* diffararK* in radiarK* without any atmotpharic corraction, 
and tha dathad linat giv* th* difftranca with atmoipharic corraction* that ignor* th* finita »i7* of th* 
fiaid*. Th* dottad linat ara tha tarn* at in T igur* 7. 
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RADIANCE (I 
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SURFACE REFLECTANCE 
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Ftgurt 9. The upwtro rtdt«nc« at tha nadir for a uniform lurfact of raflactarKt (A| 
arnt for tha thraa atmotpharic turbioitiai. 
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corrMpondlrR) radlancw for unit turbidity (T ■ 1) wbh found from th« 
curvM in Fiq. 9. Iho atmoqshoric corraction appliad hara (dashad 
linas in Piq. 8), corracta sliqhtly tha radianca for • 0.4, whila 
incraasinq tha maqnituda of the diffaranca for ■ 0.0, and havinq a 
mixad affact for Ag^ ■ 0.1. 


Iha affact of tha backqround raflactanca on tha chanqa in tha 
radianca ia shovai in Piq. 10. It is saan that tha ralativa diffaranca 
incraasas almost linaarly as tha backqround raflactanca daparts from 
tha valua of A^. Ihis linaarity can ba damonstratad with tha pravious 
aquations by first substitutinq Bqs. (2) and (3) in (1) to find tha 
r^K)ianca above a nonuni form surfacat 


I - 


■o * r 


- tA 


, -T 
(e o 


al 




1 ■ t, r. 


( 6 ) 


The ralativa diffaranca in tha radianca above tast and refaranca araas 
is qivan by: 


1 






(a“^o 


at 


t 

1 - tA 


(7) 


ar 


where tha raflactancas of tha tast and rafarance araas are the same 
(A^ ■ A^), but tliair backqrourxl raflactancas differ. Since ^A^^ « 1 
and tA « 1, tha difference is close tot 

Al 


M 

I 






‘V ■ 


( 8 ) 


As a result of tha aver^k 3 e raflactancas A^^ (for tha raqion of the 
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0 0.2 0.4 0 0.2 0.4 0 0.2 0.4 

A Bt A Bt A Bt 

Figurt 10. Tbt rdativt cbangt m radianca a* a furKtion of background raflactanca (Ag,K 
for tba thrw atmoiphartc turbrdititi. A lolid lina givat rasuUa for no atmotpharic 
cotraction, and a dathtd lina givat ratciti with uniform atmotpharic corraction. Both 
tha tatt and rafaranca araat ara of adge langth X, ■ X, ■ 0.2Q km, and A, ■ A, ■ 0.2. 
Tha background raflactivitv of tha rafararKa it Ag, "0.1. Tha dottad linat giva tha 
t6 arnf 110% tolar anca valuat of tha diffararKa in radianca. Nota that for T ■ 1, 
tha dathad lina matchat tha tolid lina. 
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t*st aquara) and (for tha rafararx^a aquara) balnq llnaarly 
dapandant cn thair background raflactanca, tha llnaar dapandanca o£ 
tha diffaranca batwaan radiancas (AI) followai 







( 9 ) 


ITm atmoapharic corraction, basad on a uniform aurfaca, la shown by 
tha dashad linaa in Pig. 10. lha atjnoapharic corraction, whila having 
no affact for T « 1, dacraaaaa tha radiancas for T ■ 0 and incraasas 
tham for T ■ 3. 


Whan tha turbidity ■ 3, tha coriactad intansitias 11a within tha 

^ 5 percant tolar anca limits only for 0.12 < < 0.18. lha 

diffaranca batwaan radiancas for tha tast and rafa''anca fialds for 

additional permutations of thair paramatars is plotted as a function 

of tha aerosol optical thickness in Pig. 11. Iha continuous lines 

show tha relative difference. Iha dashad lines result after tlia test 

radiancas are corrected for atmosperic affects that are based on the 

assumption that the surface reflectance is uniformly the same as for 

the tast field. The raflactar s of the raferervre and test fields are 

the same. Ihe reference parameters are constant (A^ ■ 0.2, Ag^ ■ 0.1 

X •0.29 km, and T -1). 
r r 

Only the backgrourvj reflectance (Ag^) decreases for each 
succeedir>g row. 0>ly the target size increases from t^ie first to the 
second column. 
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^/le ^/le ^/l« Z^/l« 


Abi - 0.4 X, - C.29 KM Ast - 0.4 X, - 1 .9B KM 



0 0.25 0.5 0 0.25 0.5 


Figure 11. The relative change in radiance as function of the aerosol optical thickness 
without (solid line} and with uniform atmospheric correction (dashed line). Each row of 
graphs is given for a different ratio of Ag,/A, which changes fo' the respective rows from 
the bottom (0, 0.6. 1.0, 2.0). The left column is for A, ■ 0.2, Ag, ■ 0.1, X, ■ 0.29 km, 
T, ■ 1, A, ■ 0.2, X, ■ 0.29 k~i. The right column is the same but for X, ■ 1.98 km. The 
dotted lines give the 15 arsd 110% tolerarKe values of the difference in radiance. 
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A (f*w i|t>r>«rAl ix>nrUuilona cam tw mAilt* r*>)Ardini| Flq. lli 


A. Hw* rAdlArw dif (•r«fH'«A Ar» ti«»iAtlvA» v4mmi th» lv»okqroirvi 
teflAi'tArK'* la totAlly blAt'k (A^^ ■ 0; bottom row). 

b. r»t lAOtAiiraa Ar«* tlvp HAm«» (Al • 0), vdu»ra th«» raf AfAnc# 

Anil taat raql^na Ara th« In th«» third riMrf (A^ ■ A^ • 0.2 Arvi Aj^^ 

■ A^^^ • 0.1). ’Hit* rAdlantw dlffi»r»m'» lApamia allqhtly on turbidity 
in th« third r»ni», whf»ri* th«* rAfaranot* Anil taat fialda Ara tin* HAim* 
alia (latt column). TTu* dltfararu'a inotaAaaa to 8 parcant v4a*ra t)w 
fiald Hizoa Ara diffarant. 

c. Hu* Aljm>.Hphar Ic oi>rraotlon (or a uniform surtAoa (dAHi)ail 
Unas) iV>a8 not v-hAnqa tha dKfaranoa iH»t\«A»an rAillancaa of tha 
rafaranoa aiuI fast AraAS, vdian t)ia turbidity T ■ 1; but tha dlffaranoa 
Im'taAsas (or T < 1 Anil dai'raAsas for T > 1. ITia corrai-taii aiuI 

itiixir rai'tail rAiiiAm'a diffaratK'aa stww tha SAm«' tarxiam'y to I la outaida 
of tha tolarAni'a limits. Ttiasa limitail liAtA IrxIlcAta t)>At this 
Atmos^)tu»r Ic corraotlon dapatulit^q ort turbidity doas not railik'a tha 
At^si^luta dlffaram'a tw»tvA>an t)wr rAillAm'as of tha tast Anil rafaram'a 
(lalds, sufflclantly to mska a siqnlfic*Ant Impiovamant in clASslf IcAtlon. 

5 . 01^‘V I UIN ArjiW 

In this saction va> simulAta tha Atjnospl^aric affact on tha 
rAdlAnoa drtactad Abova Aqrlcultural raqiona conslstlnq of tha 15 
clAssas qivun in Tabla 1. 1^a simulatad dAtA Ara clAssifiail, Atid tha 
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T«k))r li SurfAcr rrf IcK't Nnrv uiiihI In thr iiiMulAtion 



BAND 

WAVRl.RNt 

AND 

;TH (nm) 




4 

•t 

b 

7 


SURFACK 

SOO-bOO 

h(H)-7olnr 

700-800 

800-1100 

RRFRRENCES 

Soil 

0.12H 

0. 183 

0.213 

0.2 SO 

Oeerlng (1981) 

Alfalfa 30 kg/tia 

0.117 

0.11^ 

0.228 

0.30S 

Peering (1981) 

Alfalfa 4M) kg/lui 

O.OB3 

0.0H«» 

0.2^C 

0.317 

Peering (1981) 

Alfalfa lb*0 kg/ha 

O.OM 

0.044 

0. 30S 

0.4 OS 

Peering (1981) 

Alfalfa 22HO kg/lia 

O.OM 

0.033 

0.310 

0.428 

Peering (1981) 

Alfalfa 3^^0 kg/lia 

O.OM 

0.033 

0. 3b0 

0.S17 

Peering (1981) 

Alfalfa 3HM) kg/lui 

O.OM 

0.033 

0. 390 

O.SSO 

Peering (1981) 

Water 

0.010 

0.009 

0.004 

0.003 

Mar ah and l.yon (1980) 

Savanah 

0.074 

0.118 

0.142 

0.182 

Krlebel (1977) 

Bog 

0.02H 

0.038 

0.090 

0.147 

Krlebel (1977) 

I’aature land 

0.03 3 

0. 120 

0.2 SO 

0. 380 

r.rlebel (1977) 

ConlferouM land 

0.014 

0.024 

0.044 

0.11b 

Krlebel (1977) 

Corn 

0.072 

O.ObO 

0.320 

0.S4O 

Barker (1981) 

Soybean 

0. I4^ 

0.174 

0. 3h0 

0.4b0 

Barker (1981) 

Wheat Stubble 

O.Obb 

0. 100 

0.134 

0. 17S 

Barker (1981) 


Table 1-2: The normal ictnl laJlaneeN are Klven alH>ve the 

atmoMpheric In the direction of the reference areaa, 
which are uaed for c laaH If Icat Ion . The atmoapherlc 
turbidity T • 1, and the ed|;e length of 'he retennice 
field of X • 1.23 km. The background tor each 
reference area la aoll. 


SURFACK 


BANP ANP 
WAVKl.KNCTH (nm) 


4 

5 

b 

7 

SOO-bOO 

bOO-700 

700-800 

800-1100 

So 1 1 

0. Ib49 

0.20b8 

0.23b7 

0.2bS3 

Alfalfa 30 kg/lta 

0.1 SS9 

0. 182S 

0.24bS 

0.31S8 

Alfalfa 9b0 kg/lta 

0. 1284 

0. 1258 

0.2bb2 

0.32b8 

Alfalfa IbSO kg/lta 

0. 11(»8 

0.087S 

0.31S8 

0.4087 

Alfalfa 2280 kg/lui 

0. 1108 

0.0782 

0.3204 

0.4302 

Alfalfa 3bb0 kg/lta 

0. 1 108 

0.0782 

0.3bS9 

0.514b 

Alfalfa 3850 kg/ha 

0. no8 

0.0782 

0.393S 

O.S4b2 

Wat er 

0. 0702 

O.OS79 

0.0S05 

0.0438 

Savanah 

0. 1212 

0.1505 

0.1703 

0.2035 

Bog 

0.845 

0.0824 

0.1247 

0.1719 

Paature land 

0.0884 

0.1S23 

0.2bb2 

0.3853 

Cunlferoua foroat 

0.0734 

0.0705 

0.0849 

0.1441 

Corn 

0. 1 18b 

0. 1011 

0.3294 

0.53bb 

Soybean 

0.1788 

0.1990 

0.3b59 

0.4b0S 

Wheal Stubble 

0.1148 

0. 1351 

0.lb32 

0.1971 


2 « 




error of clAsaif ication is analyzed for its dependence on fields size, 
background reflectance, and atmospheric turbidity. 

Ihe surface reflectances used in the simulation are given in 
T.^)le 1 for waveler>iths corresponding to the Landsat bands {band 4i 
SOO-600 tw; 5: 600-700 mt 6j 700-8'J0 rzn; 7i 800-1100 rzn) . The 

reflectances are interpolated from data given in the cited references. 

atjm^spheric effects are derived from Pearce's (1977) Monte Carlo 
calculations of the radiance of finite fields, as they vtould be 
observeil from a satellite. Tbese raiiiances are rescaleil accordimi to 
the e«.)uations in ApfJ»nv1ix 4, to account for the ref leirtances given in 
Tal>le 1. In these simulations, the referem'e field arxl also the field 
to be classifies) ate square aixl surroundevi by a uniform background of 
different reflectance. 


An extrapolation of f^arce's (1977) radiances for 550 nm to other 

vMveletxiths is rei)uired, also. TTie w.welength dependence is expressed 

by means of the total, normal optical thickness of the atmosphere ( 

■ t ♦ t ). Itie optical thicknesses are calculated with respect to 

A K 

the optical thicknesses of particles (t ) am) mol€K*ules (t ) at 500 rim 
from the following relations: 


t” (A) - i” (550) (X/550)‘^*®^ 


(A) • 1^ (550) (A/550) 


( 10 ) 


where ^ is calculated from Pearce's data (1977) as v • 0.97. The 
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r*acallnq of tb« radiancas to account for dlffarant optical 
thicknesses is given in Appendix 4 . The reference square fields are 
taken here always with a size of 1.23 km, and for average atjnospheric 
turbidity of 1. Table 2 gives the radiances for the different 
reference fields with a backgrourxl ol sull. 


The olasslf ication of an unknovn (teat) field is performed by 
associating it with the reference field giving a minimum radiance 
difference defined byt 


r 


2 

k 


I 

4 


4 

Y. 




2 


k - J, 2, 


15 


( 11 ) 


vhere and are the radiances above the centtr cf the test and 
reference fields, respectively. The superscript k ranges through the 
complete set of integers from 1 to 15, corresponding to the 15 classes 
listed in Table 1. The turbidity is always one (T ■ 1) for the 15 
reference fields and their size is always ■ 1.23 km. 


The classification results are given in Pigs. 12a, b, and c. The 
table at the top of each figure gives details about the fields and 
their backgrounds used in the simulation. The atmospheric turbidity 
is given in the column labeled T. Colimn C indicates by the word 
"yss" the data after an atmospheric correction is applied. The 
correction is made for a surface of uniform reflectance. The fields 
are designated by a letter: R for reference field, T for test field, 

C for the test field after the atmospheric correction has been 
applied, and A, B for the one or twa fields that have the closest 
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REFLECTANCE 


Th« •ffact of turbidity aton* 


riEU) 


X 

MM 



CLASSiriEP 

or 

ON 

(Kh) 

■I 



5 

6 

7 

AS 

Coni faroua 
roraat 

Soil 

1.23 

■ 

1 

.125 

.118 

.126 

.171 

Whaat Stubbla 

Coni faroua 
roraat 

Soil 

1.23 

1 

■ 

.073 

.070 

.085 

.144 

- 

Coni faroua 
roraat 

Soil 

1.23 

■ 

Yaa 

.092 

.091 

.104 

.156 

Boq 

Boq 

Soil 

1.23 

1 

■ 

.004 

.082 

.125 

.172 

- 

Wliaat Stulibla 

Soil 

1.23 

n 

■ 

.115 

.135 

.163 

.197 

- 




Figurt 12. Examplm of atmocphaiic ttfact on dau) tication; a - Tha affact of turtMdity alona; b - Tha affact of 
lira alona; c - Combinad affact. and d - O'jtsification of watai. Each axampla givat tha radiancai of tha fialdt 
both m tabular tnd graphical form. In addition, tha raflactarKa of aach covar it givan. 
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„b» 

Th« of tiM alon* 


1 riELD 

X 

(KM) 

T 

1 RADIANO: KOK LANDSAT HANl 

CLASSIFIED 

AS 

or 

ON 

4 

s 

b 

7 

Alfalfa 

3280 

Coni fat oua 
foraat 

0.111 

1 

.097 

.066 

.297 

.411 

Alfalfa 1650 

Alfalfa 

2280 

Conlforoua 

Foraat 

1.23 

1 

B 

.066 

.309 

.423 

m 

Alfalfa 

1650 

Conlforoua 

foraat 

1.23 

1 

IQH 

.075 

.305 

.402 

m 




Figure 12 (continued) 
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REFLECTANCE 


“d" 

Clatsifkation of wator 


1 FIELD 

HI 

■1 

B 




CLASSIFIED 

w 

or 

CM 

USB 

5 

6 

7 

Watar 

Soyhaan 

0.111 

B 

■ 

.147 

.127 

.164 

.159 

Whaat Stubbla 

Watar 

Soil 

1.23 

D 


.070 

.058 

.050 

.044 

- 

Watar 

Soybaan 

0.111 

n 

Yas 

.115 

.101 

.145 

.144 

Bog 

Bog 

Soil 

1.23 

D 


.084 

.082 

.125 

.172 

- 

Whaat 

Stubbla 

Soil 

1.23 

■ 

■ 

.115 

.135 

.163 

.197 

- 



BAND 



BAND 


Figurt 12 (contimivd) 
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spBctral ctuiractarlstlcs to flalda T and C. Tha raflactancas of tha 
fialda at tha aurfaca ara shown in lowar laft fiqura. Iha 
radiancas abova tha atmosphara ara takan from tha tabla and raproducad 
in tha lowar right figura. 

Only tha atmospharic turbidity changas batwean tha rafaranca (T ■ 
1) and tast (T ■ 3) fialds in Pig. 12a. lha adga lartgth of all fiolds 
in tha simulation is X ■ 1.23 km. Iha rafaranca and tast fialds ara 
conifarous forasts surroundad by a background of soil. Because the 
radiances above the tast field are higher than those of the reference 
field, the forest is misclassifiad as wheat stubble (B) . The 
corrected data (C) is classified as bog (A). In actuality, the 
atmospheric turbidity does not usually vary so strongly between 
reference and test fields, but v^en it does, serious classification 
errors arise. 

The field size is the only parameter that changes between the 
reference (X^ ■ 1.23 km) and the test (X^ ■ 0.11 km) fields for Fig. 
12b. The fields are alfalfa of density 2280 kg/ha with backgrounds of 
coniferous forest. The test field is still classified as alfalfa, but 
of lower density (1650 kg/ha). The size difference results in a 
misclassif ication, even though the atmospheric characteristics and the 
surface reflectances are the same for both the reference and test 
regions. 

In Pig. 12c the combined effect of turbidity, b^K;kground 
reflectance, and size is examined. The test field is v^eat stubble on 
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a background of coniferous forest, for ■ 0.111 km and ■ 3. A 
coRfMrison of the radiances of this test field with a reference field 
of k^at stubble shows that t!.e atmospheric effect is weak in this 
case due to the relatively dark background. Ihsrefore ':he field is 
classified correctly as vAieat stubble. Ihe atmospheric effect would 
be much stronger for a uniform field of \4)eat stubble; thus, the 
resultant uniform atmospheric correction to the radiances in (T) 
results in a decrease of the radiances to the result shown in (C) . 

This decrease is caused by the uniform atmospheric correction, which 
does not account for the adjacency effect. The res^ilt of this 
correction is that the field is misclassif led m bog (A). This 
example shows how uniform atmospheric correction can cause a 
misclassif ication of data that were classified correctly before. 

An example of classification of water basins is shown Pig. 12d. 

A square water pond of 0.111 km on a side with a background of 
soybeans for T ■ 3 is misclassif led as wheat stubble due to the large 
increase of radiances from line R to line T. The uniform atmospheric 
correction causes here only a slight correction, and the water pond is 
still misclassif led, this time as bog. This example shows that small 
reservoirs of water on bright backgrounds can be misclassif ied from 
space, and a uniform atmospheric correction carviot provide an adequate 
correction. 

6. CONCLUSIONS 


The simulation study presented in this paper shows that 
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nonuiitonn r*fl«ct«nc« of th« oarth's aurfaca rasults in an Important 
modification of affact of tha atmosphara on raa»ta sansing of tha 
aurfaca and thua on claaaif ication of lurfaca faaturaa. Iha adjacancy 
affact, which incraaaaa tha radianca datactad abova dark fialda 
aurroundad by adjacant bright araaa (and vica varaa) craataa an 
atmoapharic affact, which cannot ba corractad by mathoda utilizing a 
uniform aurfaca. Tharafora, futura davalopmant of atmoapharic 
corraction algor ithma haa to account for tha nonuniform natura of tha 
aarth'a aurfaca. 
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APPEiriX 1 


List of notations frsqusntly ossd in ths papsr. 


A. 


B. 


Surf SCO rsflsctancss 

A - Avsrags of rsflsctancs of both a fisld and its 

background for calculation of ths multipla 
acattarinq of light batwaan tha surfaca and tha 
atmosphara. 

Ag - Background raflactanca (raflactanca of tha araa 

surrounding a fiald). 

Ajj^, Ag^ - Background raflactanca for tha tast and 
rafaranca fialds, raspactivaly. 

Aq - Weightad avaraga of tha surfaca raflactanca of 

tha antira surfaca for calculation of tha 
diffusa trimsmission through tha atmosphara. 

A^ - Raflactanca of a given fiald. 

A^ - Raflactanca of a rafaranca fiald. 

Ag - Raflactanca of tha surface within IPOV. 

A^ - Raflactanca of a tast fiald. 


Radiances 

I* - Absolute radiance. 

I - Normalized radiance equivalent to tha apparent 

raflactanca and given by I - itIV(P F ) . 

o o 

- Radiance transferred diffusively through the 
atmosphere from tha surfaca to a datactor. 

I • Ihe radiance which would be measured if no atmosphere 

o 
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v«r« pr«Mnt 


Iq - R«diai>c« of light acattotod tron tha direct aolar 
ba«n Into tha datactor fiald of viaw v.ltiK>ut balng 
acattarad by tha aurfaca. 

Ig - Radlanca of light raflactad by tha obaarvad fiald and 
diractly tranantttad through tha atmoaphara. 

- Radianca abova uniform aurfaca. 

• Radianca abova nonuniform aurfaca. 

C. Othar notationa 

- Atmoapharic turbidity abova tha rafaranca fiald. 

- Atmoapharic turbidity eibova tha taat fiald. 

- Edga langth of tha rafaranca fiald. 

- Edga length of tha taat fiald. 

- The total transmiasion of tha sunlight to tha 
o 

aurfaca. 

D • Diffusa transmission through tha atmosphere. 

A M 

- Total atmospheric optical thicKnass (i t ) . 

- Aerosol optical thic^ nass. 

- Mr.lacular optical thickness. 
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APPENPIX 2 - APPLICATION OP DAVE'S CODE 


Th« equations giving the radiance above the center of an 
infinitesimally small square lying on a uniform background are 
derived. Chandrasekhar (1%0), gave an equation for the upward 
radiance above an atmosphere on a uniform surface that reflects light 
according to Lambert's law: 


I 


u 



(EA •r DA) 


( 2 - 1 ) 


where E is the direct upward transmittance by the atmosphere. In the 
case of a small object of reflectance, A^ laying on a background of 
reflectance Ag, the upward radiance is given by: 


I (A^, Ag, T, X -► 0) 



(EA^ DAg) 


( 2 - 2 ) 


Equation (2-2) is obtained from Bq. (2-1) by substituting A^ for the 
reflectance in the term of direct transparency, and Ag for terms of 
diffuse transparency and reflection. Ihis can be done as long as the 
object is not big enough to affect the diffuse component (edge length « 
300 m) . Dave's results for an infinitesimal field and an infinitesimal 
field of view of the detector are compared with the Monte Carlo 
calculations for a finite field of view of footprint of 30 m and object 
field of edge length of 40 m (Figs. 5 and 6). The differen^'es shown in 
Figs. 5 and 6 are partially due to the finite sizes of the footprint and 
the object field and in addition can be due to computational errors or 
small differences between the atmospheric model used in the Monte Carlo 
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calculationa and In the Dave calculatlona 


The imall differences ilo not Introduce significant errors in the 
adjacency effects^ because of the wide range of latter (Mekler ami 
Kaufman, 1980). *n>us, Pearce's data are reliable for studies of the 
adjacency effects. 
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APPENDIX 3 - PHOPORTIGNAL RADIANCE TRANSFX3RMATION 


ITie •quAtlona thAt ar« lUMd to transform th« radiances from the 
values in Piqs. 5 and 6 for one set of surface reflectances to new 
radiaiKres for different sets of surface reflectances are derived in 
this appendix. ‘Ihese equations are use to calculate the data on Piqs. 
7-8 and 10-11. The radiance at the center of a square of reflectance 
A^ of any edqe lenqth (X) can be formally qiven by (see Bq. 5): 

Y 

I (Aj. T, «)-!„* r- T A ~ '“t * “n’ (3-‘> 

A 


where E > (exp (-t )), A and A,, are reflectances obtained as wciqhted 

O A w 

average of the reflectance of the square and its backq>‘ound: 


A 

a 


rt ♦ 

a t 


(1 - a ) 

A 


■ a. 


a - 


a ) 

a 


(3-2) 


where a and a depend on the atmospheric optical characteristics. 

a u a 

and can be expressed as complex solutions to the r^Kliative transfer 
equation, we shall not find these solutions here but express by 
using the Monte Carlo data; and shall roughly estimate the value of 
a . Rough estimation of a does not cause a significant error in Uq. 
(3-1), since the term (?A ) is relatively small. In the present work 
0.2, affecting the final radiance by 0-10 percent. Thus, an error 
in a of 10-20 percent causesOan error in the radiance of 0-2 percent. 

A 
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trarvtfomwitlon from known valuas of I for a sot of 
roflectancoa A^) to new values of I for another set (A^'» A^') is 
qlven hyi 


(a|, X) - (1 j^(A^, K^, T, X) - ij • R^) 


1 - *A 

■- lA R ’ ^ *0 

a a 


(3-3) 


The (lata In Figs. 7-8 and 10-11 are calculated for: 


R 


A 



I 



% 


(3-4) 


This racio has been used to obtain the simplified expression (3-3). 
In this transformation, the linear dependence of the radiance on 
surface reflectance is transformed exactly, v4iile on the nonlinear 
part (1 - iA ) an approximate transformation was applied using the 
following expression for A.: 

A 



(3-5) 


wt^ere X is the edge length. The value of R ■ 6 km is found to give 
good results and is used in the calculations. This approximate 
expressions might introduce errors in A of (AA ■ ♦ 0.1). Since the 
value of i used in the present paper is t £ 0.21, the resultant error 
in radiance is less than 2 percent. 
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APPE343IX 4 - GENERAL RADIANCE TRANSPGRMATJQN 


In this appendix we shall explain the transformation applied to 

the radiance I(A^» Ag, x) above a square of reflectance A^ and 

edge length X lying on a background of reflectance Ag and for optical 
A M 

thicknesses t and t to obtain another r^Kiiance for a new set of 
parameters A^*, Ag*, and giving I(A^*, \g*, X). 

ITie transformation is done in two steps. First, the transformation is 
nude with respect to surface reflectances and then transformation for 
the optical thickness. Ihe transformation for the surface 
reflectances is done by Eqs. (3-1) and (3-2). Substituting Ag from 
Bq. (3-2) into Bq. (3-1), we find that : 


Y 

I (Aj, Ag, t", X) . A^ (E ♦ DOp) ♦ 

-«d' 


(4-1) 


For a uniform surface of reflectance A^, we can rewrite Eq. (2-1) in 
the form: 


(A^) 


1 - tA. 


(EA^ 


DA^) 


(4-2) 


From Eqs. (4-1) and (4-2) we can write: 

I (A^, Ag, t”, X) - l“ (A^) + YjjE + C ~ 


(4-3) 


v4iere: C ■ yd (1 - a ) 

O D 
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I 

■ (1 - tA^) (1 - lA^) 


, I - A ) 

1 ♦ -4 2 _ 

1 - tA^ 


Expression (4-3) can be sinnplifled by neqlectinq terms of the order of 
(*^Aa), TTiis approximation was found to introduce an error of 0-4 
percent in the resulting radiances for the present calculations. Ihe 
resultant coefficients are G ■ 0 and G 2 " 1. Hence 


I (A^, 1^, t”, X) 


1“ (A^) ♦ c 



(4-4) 


'nie transformed radiance: 


I (A^*. A^*. 



(A^) ♦ C 



(4-5) 


is calculated by using C from Bq. (4-4) tor the known value of I(A^# 
A«, X), and using A and A^' from Bq. (3-5). The value of 

* WB3 calculated from Iq and know value of by Bq. (4-2). 

The transformation on the optical thickness is based on the 
following assunptions. 


a. ITie radiance is both linear on the aerosol and molecular 
optical thickness. 
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b. It Is possible to separate the dependence of the radiance on 


the aerosol optical thickness from its dependence on the molecular 
optical thickness. 


c. Mathematically, these assijnptions allows us to express the 
the rexliance byt 

I t”) ■ io ”*• 


Table (4-1) shows examples of the radiances calculated by Bq. (4-6) 
from radiances at other optical thicknesses, and compared with 
radiances obtained directly. Ihe edge length of the square (X) for 
v^ich the reidiance is given is also tabulated, in addition to the 
values of an Ag. 


Equation (4-6) can be used to transform the radiance to new 
values of optical thickness by: 


I (t*‘, ,"•) 




M* 


I (0, T ) - A T_ 


(4-7) 


I (T^, t”) 


M. 


- I (0, T ) t 


~ + A. 




I 



T«bl« 4-l i Comparison between rsdisnees obtained by Plqs. (4-6) 
and radiances calculated directly. 



RADIANCE! 

i USED IN 

riGS.(4-6) 


T^0 


T**-0.10 

T^-0.21 

RESULTANT 

RADIANCE 

t”-0.1. T^-0.637 

RADIANCES 

OBTAINED 

DIRECTLY 

X(KM) 

> 

o 


T^0 

1.23 

0.4 

0.2 

0.4 

0.391 

0.381 



0.361 

0.366 

0.179 

0.6 

0.0 

0.6 

0.544 

0.479 

0.349 

0.354 

0.47 

0.2 

0.4 

0.2 

0.229 

0.256 

0.310 

0.315 

0.758 

0.0 

0.6 

0.0 

0.063 

0.120 

0.234 

0.250 

oo 

0.0 

- 

0.0 

0.039 

0.051 

0.075 

0.088 

oo 

0.2 

- 

0.2 

0.218 

0.231 

0.257 

0.251 

oo 

0.4 

- 

0.4 

0.405 

0.408 

0.414 

0.420 

oe 

0.6 

- 

0.6 

0.596 

0.601 

0.611 

0.605 


45 










REFERENCES 


Barker# J, 1981: Private communication. 

Chandrasekhar# S.# 1960: Radiative transfer. Dover# New York# p. 

273. 

Dave# J. V. and J. Gazdag# 1970: A Modified Fourier Trarwform Method 

for Multiple Scattering Calculations in a Plane-Parallel Mie 
Atmosphere# Appl. Opt. # 9# 1457-1466. 

Dave# J. V.# 1980: Effect of atmospheric conditions on remote sensing 

of a surface nonhomoqenity. Photoq. Ehg. and Dam. Sen.# 46# 1173. 

Fraser# R. S.# 1974: Computed atmospheric corrections for satellite 

data. Proc. of the SPIE Conference on Scanners and Imagery 
Systems for Earth observation# SPIE# ^# 64. 

Fraser# R. S.# Bahethi# Om P. # and A1 -Abbas# A. H. 1977: The effect 

of the atmosphere on classification of satellite observations to 
identify surface features. Remote Sensing of Environment# 6# 229. 

Griggs# M.# 1973: A method to measure the atmospheric aerosol content 

using ERTS-1 data. Proc. of the Third ERTS-1 symposiim# NASA 
Sp-351# pp. 1505. 

Griggs# M# 1979: Satellite observations of atmospheric aerosols 

during the COMET cruise. J. Atmos. Sci.# 36# 695. 

Haba# Y.# Kawata# Y.# Kusaka# T.# and Uneo# S.# 1979: The system of 

correcting remotely sensed earth Imagery for atmospheric effects. 
Proc. of the 13 Int Symp. on Remote Sensing of Environment# Ann 
Arbor# Michigan# pp. 1883. 

Herman# B. M. and Browning# S. R.# 1975: The effect of aerosols on 

the earth-atmosphere albedo. J. Atmos. Sci.# 3i2# 1430. 

Horvath# R. B.# Polcyn# J. G.# Fabian# C. # 1970: Effects of 

atmospheric path on airborne multirpectral sensors. Remote 
Sensing of Environment# U 203. 

Kaufman# Y. J. and Joseph# J. H.# 1981: Evaluation of surface albedos 

and extinction characteristics of the atmosphere from satellite 
iamges. Accepted for publication in J. Geophys. Res. 

Kawata# Y.# Haba# Y.# Kusaka# T.# Terashita# Y.# and Ueno# S.# 1978: 
Atmospl^re effects and their correction in airborne sensor and 
Landsat MSS data. Paper presented at Proceedings of the 12th 
International Symposiun on Remote Sensing of Environment# ERIM# 

Ann Arbor# Michlg^^n. 

Kriebcl# K. T.# 1977: Reflection properties of vegetated surfaces: 


46 



tables of measured spectral biconical reflectance factors, 
Uhiversitat Nunchen, Wlssenschaftliche Mltteilung, No. 29, 

Munchen, W. Germany. 

Labovltz, M. L., Toll, 0. L., and Kennard, R. E., 1980t Preliminary 
evidence for the influence of Physiography and scale upon the 
autocorrelation furtction of remotely sensed data, N^SA Ifl 82064, 
NASA/Goddard Space Plight Center. 

Marsh, S. E. and Lyon, R. J. P. , 1981 t Quantitative relationships of 
near surface spectra to Landsat radiometric data. Remote Sensing 
« of Environment, j^, 241. 

Mekler, Y. and Kaufman, Y. J., 1980: Ihe effect of earth's atmosphere 

on contrast reduction for a nonuni form surface albedo and two 
halves field. J. Geophys. Res., 4067. 

Mekler, Y., Quenzel, H., Qhring, G., and Marcus, I., 1977: Relative 

atmospheric aerosol content from ERTS observations. Geophys. 

Res., 967. 

Otterman, J. and Fraser, R. S., 1976: Earth-atmosphere system and 

surface reflectivities in arid regions from Landsat MSS data. 
Remote Sensing of Environment, 247. 

Otterman, J. and Fraser, R. S., 1979: Adjacency effects on imaging by 

surface reflection and atmospheric scattering: Cross-ir radiance 

to zenith. Appl. Opt., 2852. 

Pearce, W. A., 1977: A study of the effects of the atmosphere on 

thematic mapper observations. Rep. 004-77, EXl&GAiashington Anal. 
Serv. Center, Appl. Sys. Dep., Riverdale, Maryland, 136 pp. 

Peterson, J. T., Flowers, E. C., Berri, G. J., Reynolds, C. L. , and 
Rudisill, J. H., 1981: Atmospheric turbidity over Central North 

Carolina. J. Appl. Met., 20, 229. 

Pitts, D. E. and BeKihwar, G., 1980: Field size, length, and width 

distributions based on LACIE Ground Truth data. Remote Sensing of 
’ Environment, 201. 

Potter, J. F., 1976: Correcting leuidsat data for changes in sun 

angle, haze level, and background reflectance. Proc. of the 1976 
machine processing of remotely sensed data symposiun. IEEE 76CH 
1103-IMPRSD, 2B-6. 

Richard, P. R., 1979: Determination of noise equivalent reflectance 

for a multispoctral scanner - a scanner sensitivity study, NASA 
Tech, paper, NASA Lyndon B. Johnson Space Center, 1575. 

Richard, P. R., 1981: Private communication. 

Schowengerdt, R. A. and Slater, P. N., 1979: Literature survey and 

review of atmospheric effects in remote sensing and their 


47 


influenc on clM.lflc.tlon. OBI, Inc. 1400 Sprln, StrMt, Sllvt 
Spring, Maryland, 20910. 


Spring, MarylaivJ, 20910, 

Symnposiun on Remote Sensing of Eiivlronmenc, ERIM, October, » 
pp. 895-934. 

TMrner. R. E., 1978t Elimination of atmospt.erlc effects from rem^e 
sensor data. Proc. of the 12th International Sywposliw on Remote 
Sensing of Eiivironment, ERIM, pp. 783. 


9 


48 


